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Abstract 

A subset of neutrino- induced upward through-going muons in the Super-Kamiokande 
detector consists of high energy muons which lose energy through radiative processes 
such as bremsstrahlung, e + e~ pair production and photonuclear interactions. These 
"upward showering muons" comprise an event sample whose mean parent neutrino 
energy is approximately 1 TeV. We show that the zenith angle distribution of upward 
showering muons is consistent with negligible distortion due to neutrino oscillations, 
as expected of such a high-energy neutrino sample. We present astronomical searches 
using these high energy events, such as those from WIMP annihilations in the Sun, 
Earth and Galactic Center, some suspected point sources, as well as searches for 
diffuse flux from the interstellar medium. 
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1 INTRODUCTION 



In order to select neutrino events with the highest energies, we consider muon 
neutrino interactions in the rock around the detector, because the effective 
target volume is very much increased [1]. To separate neutrino-induced muons 
from cosmic ray muons, we select only upward-going muons, since the back- 
ground from downward going cosmic ray muons overwhelms any neutrino- 
induced muons from above. Muons penetrating the detector have energies of 
at least several GeV and point along the neutrino direction within a few de- 
grees, allowing astrophysical studies. Neutrinos originating from cosmic point 
sources are expected to have harder energy spectra than the background of 
atmospheric neutrinos. Some of these high energy neutrino-induced muons un- 
dergo radiative energy loss. We identify these muons as "showering muons". 
High energy muons are correlated with high energy neutrinos, which allow us 
to statistically select an extremely high energy parent neutrino sample from 
the Super-Kamiokande data. Our aim is to extract this sample for physics and 
astronomy studies. 



The Super-Kamiokande (Super-K) experiment [2] has previously analyzed 
two topologically distinct categories of upward muons caused by neutrino 
interactions in the rock below the detector: muons which exit the detector 
(called "through-going") and those which stop inside the detector (called 
"stopping" ) [3] . The parent neutrino energy of upward stopping and through- 
going muons for atmospheric neutrinos is peaked at ~ 10 GeV and ~ 100 
GeV respectively. In 1646 days of data, Super-Kamiokande detected 1856 up- 
ward through-going and 458 upward stopping muons. Oscillation results using 
only upward through-going muons and stopping muons have been reported in 
Refs. [3] and [4] . Both upward muon samples have been combined along with 
the contained events to do the most precise atmospheric neutrino oscillation 
studies [5,6]. Moreover, astrophysical searches for annihilation signatures of 
WIMPs [7], point and diffuse sources [8], and neutrinos from GRBs [9] have 
been done with upward muons. 



In this paper we bifurcate the upward through-going muon sample into upward 
showering and upward non-showering events. This classification is done on 
the basis of the physical mechanism for muon energy loss above 1 TeV. The 
showering sample is then analyzed to see what these highest-energy events 
might reveal. 
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Fig. 1. Average muon energy loss/tracklength for about 100 simulated muons at 
various energies from 10 GeV to 30 TeV. The theoretical curves for energy loss due 
to ionization and other radiative processes have been obtained from Ref. [11]. 



2 MUON ENERGY LOSS 



At low energies, from 100 MeV to 100 GeV, muon energy loss is dominated by 
atomic ionization as described by Bethe-Bloch; the minimum ionization en- 
ergy loss in water is 1.99 MeV/cm [10]. At higher energies, radiative processes 
such as e + e~ pair production, bremsstrahlung, and photonuclear interactions 
become important, resulting in catastrophic energy losses with large fluctu- 
ations. In water, the critical energy where ionization energy loss is equal to 
the average radiative energy loss is 1.03 TeV [11]. The relevant processes have 
recently been reviewed in Refs. [12,13]. 

To develop an algorithm for showering/non-showering separation and test its 
efficiency, we used a GEANT 3 based detector simulation to simulate muons at 
several energies between 10 GeV and 30 TeV. We checked this Monte Carlo 
sample to make sure that the muon energy loss in our detector simulation 
agrees with calculated valuesfll]. The average energy loss/tracklength at each 
muon energy is shown in Fig. 1, in good agreement with the calculation. 

All details of the muon energy loss calculations for different processes in 
GEANT 3 are discussed in Ref. [14]. The cross-section for bremsstrahlung is 
obtained from Ref. [15]. The energy loss for direct electron-positron pair pro- 
duction is obtained from Ref. [16]. Cross-sections for photonuclear interactions 
are obtained from Ref. [17]. 
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3 DATA AND MONTE CARLO USED 



For this paper, we used all upward through-going muons from the full Super- 
K-I dataset from May 1996 to July 2001 spanning 1645.9 days of data. We 
require that reconstructed tracks have a minimum length of 7 m, which re- 
duces the contamination from photo-produced upward pions to 0.01% [18,19]. 
This path-length cut of 7 m corresponds to a muon energy threshold of 1.6 
GeV. As discussed in Ref. [5], all upward muon events require between 8,000 
and 1,750,000 photoelectrons. A specialized analysis looking for upward muons 
in events with greater than 1,750,000 photoelectrons has been reported else- 
where [20]. 

We applied the showering algorithm to simulated upward through-going muons 
from the standard Super-Kamiokande atmospheric neutrino Monte Carlo sam- 
ple [5] . The input atmospheric neutrino flux used for this Monte Carlo sample 
was from Ref. [21] for neutrino energies up to 1 TeV. At 1 TeV, the calculated 
flux in Ref. [22] was rescaled to that in Ref. [21], and used for neutrino en- 
ergies up to 100 TeV. The uncertainty in the absolute atmospheric neutrino 
flux is about 20% below 1 TeV and and ranges from 25% [23] to 40% [24,25] 
around 1 TeV. The atmospheric neutrino flux in Ref. [23] is about 15-20% 
greater at 1 TeV as compared to Ref. [21]. Neutrino interactions with the rock 
outside the detector (assumed to be "standard rock") as well as in the water 
were simulated using the NEUT [26] interaction package. Neutrino interactions 
in the rock were simulated up to 4 km from the center of the detector. The 
total livetime simulated was equal to 100 years. 



4 LIGHT PRODUCTION MODEL 



Below the critical energy, an ionizing muon produces a constant amount of 
Cherenkov light per unit tracklength. However, a muon which undergoes radia- 
tive energy loss produces high-energy photons which create electron-positron 
pairs, thus increasing the total Cherenkov light in the detector. If we can 
calibrate the total Cherenkov light emitted by a normally ionizing muon (af- 
ter accounting for the various sources of light attenuation in Super-K) then 
any electromagnetic shower associated with the muon will emit excess light 
over this amount. Thus we need to calibrate the total "luminosity" of ionizing 
muons with energy well below the critical energy of 1 TeV. 

Using the muon entry point and direction as inputs, and accounting for the ef- 
fective water attenuation length (L att ) and geometrical acceptance corrections, 
we apply various corrections to the raw PMT photoelectrons. The corrected 
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number of photo-electrons of each PMT in the Cherenkov cone is: 



( Jul.) 

g corr (corr. pe) = K , (1) 

where q raw is the raw number of photoelectrons detected by each PMT; d w 
is the distance traveled by the photons from the point along the muon track 
where the photon is emitted to the PMT which detects it; F(9) accounts for the 
PMT angular acceptance and shadowing; and K is a normalization constant 
(= 1/2500 cur 1 ) which makes the corrected photoelectrons the same order of 
magnitude as the original raw photoelectrons. With these corrections to the 
charge of each PMT, the units of q corr and other terms obtained from q corr are 
corrected photoelectrons. 

We then calculate the average charge in small tracklength intervals (50 cm) 
along the muon track: 

N pmt 

Qcorr ^ ^ Qcorr j (2) 

J V* k=i 

where N pmt is the number of PMTs whose projected distance along the muon 
track is within the 50 cm path-length interval. The statistical uncertainty in 

Qcorr is S iven b y : 

N-nrat I _ \ 2 



1 JVpmi / 

2 L \Hcorr 

corr /V^ ^ * n 



(3) 



iraw 



assuming that the error in raw PMT charge arises from the square-root of 
the number of photoelectrons. We then plot Q % corr along the muon track. An 
example plot for a typical ionizing Monte Carlo muon event can be seen in 
the top panel of Fig. 2. We also calculate the average of Q % corr along the muon 
track as follows: 

N-3 , N 
E [Qiorr/VQUJ 
{Qcorr) — JfZ^ , (4) 

E i/rii 

■ ^ corr 

where indicates the number of 50 cm bins in which the muon track is sub- 
divided. In Eqn. 4, the sum goes from 4 to iV — 3, since the first and last 
1.5 m from the muon track are excluded. This is because corrections for PMT 
acceptance and shadowing are not accurately modeled near the wall. 

In addition, because of the effects of light scattering into the Cherenkov cone, 
longer path-length muons have greater corrected mean charge compared to 
short path-length muons. This is accounted for by using a path-length correc- 
tion to the average charge of an ionizing muon. The corrected charge distribu- 
tion of a typical ionizing muon is compared to that of a showering muon with 
the same entry point and direction in the bottom panel of Fig. 2. Our aim is to 
construct a function which quantitatively distinguishes between these based 
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Fig. 2. The corrected photoelectrons for an example simulated ionizing muon with 
energy 20 GeV (top) and a showering muon with energy 10 TeV (bottom) simulated 
with the same entry point and direction. 

on shape and total collected photoelectrons. 



5 SELECTION ALGORITHM FOR SHOWERING MUONS 



We define a showering x as follows: 

I \y£ cor r 
i=4 I °Q 



2 _ \ ST^ j [Qcorr (Qcorr)] I 

Xshawering "^.gL ™ ^ f ' K ) 



where N is same as in Eqn. 4; (Q corr ) is defined in Eqn. 4; <JQi corr is the 
statistical error in Q l corr and defined in Eqn. 3. Equation 5 measures deviations 
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in the histogram from a flat distribution. 



About 0.5% of events in the Monte Carlo consist of through-going muons with 
bad fits or stopping muons which are misclassified as through-going muons. 
These events show large values of \ 2 ey en though the total number of pho- 
toelectrons in such events is very small. To eliminate these cases, we define 
another variable to distinguish between a showering and non-showering muon: 



where (Q CO rr) is defined in Eqn. 4 and Q eX p(l) is a path-length dependent es- 
timate of the expected charge of a normally ionizing muon. The magnitude of 
A indicates how much more Cherenkov light is present relative to a normally 
ionizing muon. The comparison of these variables for data and 100 year at- 
mospheric neutrino Monte Carlo is shown in Fig. 3. The distributions of these 
variables for data and Monte Carlo is in good agreement. 

We applied this algorithm to the 37287 upward through-going muon events in 
the 100 year neutrino Monte Carlo [5]. The scatter plots of the distribution 
°f xlhowering an d A are shown in Fig. 4 for a 10 year subset of the 100 year 
Monte Carlo. They are presented separately for samples of "true" showering 
and non-showering muons, which are defined based on the true muon energy 
loss per unit tracklength AE/L (where AE/L is the true muon energy loss 
in the inner detector for a muon of path- length L), with the separation value 
equal to 2.85 MeV/cm. 

The separation between showering and non-showering events is decided by 
the line shown in Figure. 4. With this cut we found a total of 5747 upward 
showering events in the Monte Carlo. We then calculate the efficiency, defin- 
ing true showering muons as those with AE/L > 2.85 MeV/cm. This cut 
selects about 70% of such events, while the misidentification of events with 
true AE/L < 2.85 MeV/cm is about 5%. These cuts provide reasonable purity 
and background contamination in forming a sample for astrophysical studies 
and neutrino oscillation. 

The parent neutrino energy spectra of all upward stopping, non-showering 
through-going, and showering through-going muons are shown in Fig. 5. The 
mean parent neutrino energy of upward showering muon events is peaked at 1 
TeV. Thus the showering dataset constitutes the highest energy neutrinos seen 
in Super-K. The total number of upward showering muons, is approximately 
one-fifth the total number of upward through-going muons. 

The estimated muon angular resolution of upward showering muon events is 
about 1.4°, where angular resolution is defined as the average angular separa- 
tion between true muon direction and reconstructed muon direction. The cor- 
responding angular resolution for upward stopping and through-going muons 
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Fig. 3. Comparison of showering variables as described in Eqn. 5 and 6 for Monte 
Carlo and 1645.9 days upward through-going muon data. The solid line is normalized 
to total number of events. 

is about 2.4° and 1.3° respectively. For higher energy muons, the additional 
photons generated through radiative processes make it harder to reconstruct 
the muon track direction. Hence the angular resolution is slightly worse for 
showering muons compared to the non-showering muons. 



6 BACKGROUND SUBTRACTION 



When we applied the showering muon algorithm to the upward through-going 
muon sample from the 1646 days of data, we found a total of 318 showering 
events. The upward muon data sample is contaminated near the horizon from 
downward going cosmic ray muons which appear as upward going because of 
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Fig. 4. Distribution of A vs Xshawering f° r even ts from a 10 year Monte Carlo sample 
(which is a subset of the 100 year atmospheric neutrino Monte Carlo) with true 
AE/L > 2.85 MeV/cm (left panel), i.e. "true showering events," and true AE/L < 
2.85 MeV/cm (right panel), i.e. "true non-showering events." The solid line indicates 
the cut used to separate showering from non-showering muons. 
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Fig. 5. Neutrino energy spectra of all 3 categories of upward muons in Su- 
per-Kamiokande in 100 yrs of atmospheric neutrino Monte Carlo. 

multiple Coulomb scattering and the finite angular resolution of the fitters. 
To estimate this background in the showering muon sample, we use the same 
procedure as that applied for stopping and through-going muons [5]. We ap- 
ply the showering algorithm and select a subsample of near horizon downward 
through-going muons with < cos0 < 0.08, where O is the zenith angle 
denned with the vertical upward direction at cos© = —1. We then estimate 
the background from the downward muon sample by extrapolation below the 
horizon [5]. The zenith and azimuthal distribution of near-horizontal shower- 
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Fig. 6. Zenith and azimuthal distribution of showering muons in the region of 
— 1.0 < cosQ < 0.08. Regions (1) and (2) in bottom panel refer to thick and thin 
parts of the mountain respectively. 

ing muon events is shown in Fig. 6. The expected background from showering 
muons is estimated to be 9.15±5 3 3 °, all in the bin from < cos < —0.1. This 
background is subtracted from the upward showering muon dataset for oscilla- 
tion studies. For astrophysical studies, we assign a weight to this zenith angle 
bin which is equal to the ratio of background-subtracted to non-background 
subtracted events. 



7 OSCILLATION ANALYSIS 



Previous analyses have shown that the Super-K atmospheric neutrino dataset 
is consistent with neutrino oscillations with Am 2 ~ 0.0025 eV 2 and sin 2 29 = 



12 



100 

90 

80 

| 70 
CD 

> 60 

O 50 
CD 

-Q 40 

^ 30 
20 
10 

° -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 
COS© 

Fig. 7. The expected zenith angle distribution of upward showering events with and 
without neutrino oscillations. Both histograms are normalized to the livetime of the 
dataset. 

1.0 [5,27]. Given these neutrino oscillation parameters, the oscillation proba- 
bility is negligible for neutrinos with an energy of 1 TeV at all path-lengths 
through the Earth. The zenith angle distributions of upward showering muons 
without oscillation and with oscillated parameters obtained from Ref. [5] arc 
shown in Fig. 7. As expected, the difference between these distributions is neg- 
ligible. We will demonstrate this by comparing the zenith angle distribution 
of data and Monte Carlo. 

The chi-square function used for comparison of data and Monte Carlo is the 
same as that used for our published oscillation analysis [5], which is based on 
the pull method [28]: 




X 



10 

E 

i=l L 



N obs _ N exp 
„stat 



2 N a . 



1 2 



+ E 

3=2 



e 3 

sys 

3 - 



N s . 



nt = Nrn + E fa: 

3=1 



(7) 
(8) 



In Eqn. 7, N° is the number of observed events in i bin, N^ xp is the ex- 
pected number of events considering both oscillation and systematic uncer- 
tainties, a % stat combines statistical uncertainties in the data and Monte Carlo 
simulation, N sys is the number of systematic errors used for the fit, and N° sc is 
the expected number of events in this bin assuming to v T oscillation with- 
out considering the contribution from systematic uncertainties. The factor /j 
represents the fractional change in the predicted event rate in the i th bin due 
to a variation of the parameter €j. The second sum in the x 2 definition collects 



13 



Table 1 



Summary of systematic errors and the best fit at null oscillation used for upward 
showering muons. 



No 


Systematic uncertainty 


a{%) 


Best-fit (%) 


1 


Absolute Normalization 


Free 


20.4 


2 


K/ir ratio 


20 


-6.75 


3 


Axial vector Mass (M A ) 


10 


0.15 


4 


Multi-pion production (model-dependence) 


1 


-0.44 



the contributions from the systematic uncertainties in the expected neutrino 
rates. 



As we are doing oscillation analysis with only one data sample, and since the 
absolute normalization is considered to be free as done in Ref. [5], we only 
need to consider those systematic error coefficients for which the response is 
non-uniform in different zenith angle bins. The effect of all systematic terms 
for which // is same for all i can be incorporated in the absolute normalization. 
Of all the systematic terms used in Ref. [5] , there are only three for which the 
response is different in various zenith angle bins. These are shown in Table 1. 

The atmospheric neutrino Monte Carlo does not include the attenuation of 
neutrino flux in the Earth. We did an estimate of the effect of neutrino absorp- 
tion for the showering muons using the models in Ref. [29] and the expected 
decrease in number of events is about 1.1% . Therefore we do not incorporate 
it into the zenith angle prediction. Hence N sys = 4 in Eqns. 7 and 8. 

We vary these four ej in order to minimize x 2 f° r each choice of oscillation 
parameters sin 2 29 and Am 2 . Among these, only three contribute to the x 2 
because there are no constraints on the absolute normalization. The mini- 
mum x 2 value, Xmin = 3.54/7 DOF, is located at (sin 2 29 = 1.0, Am 2 = 
1.05 x 10 -2 eV 2 ). The number of DOF is found by 10 terms in the x 2 sum P 1US 
three systematic constraints in the x 2 sum minus four minimized parameters 
minus the two physics parameters of sin 2 29 and Am 2 . For null oscillation, 
(sin 2 29 = 0), we found a x 2 value of 6.21 for 9 DOF, where only the overall 
normalization is a free parameter. The best fit values for the systematic un- 
certainties for null oscillations are indicated in Table 1. In Eqn. 8, the best-fit 
value of e for absolute normalization for the case of null oscillation is about 
20.4%. The minimum x 2 is not significantly different than the x 2 obtained for 
null oscillation. 

Thus, we satisfy an important consistency check of the neutrino oscillation 
parameters determined in Ref. [5], namely, that this dataset (with mean neu- 
trino energy of ~ 1 TeV) is consistent with null oscillation. In future studies, 
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Fig. 8. Observed number of upward showering muons as a function of zenith angle. 
The dots represent the observed events with statistical error bars. The solid line 
shows the expectation at null oscillation evaluated using Eqn. 8 and from best fit 
values in Table 1. 

this dataset could serve as an additional high energy bin, constraining high 
Am 2 solutions as well as certain non-standard oscillation scenarios [30,31]. 
The zenith angle distributions of upward showering muon data along with the 
expected distribution from null oscillations and using the best fit values from 
Table 1 are shown in Fig. 8. 



8 ASTROPHYSICAL SEARCHES 



The high energy nature of this subset of the Super-K neutrino data provides 
an advantage when searching for astrophysical neutrino point sources, neu- 
trinos from WIMP annihilation, and diffuse neutrino flux from the galactic 
plane. The main advantage in using this high energy sample for astronomy is 
that the atmospheric neutrino background is very much reduced due to the 
steeply falling atmospheric neutrino spectrum [32,33,34,35]. Furthermore, the 
high momentum of the incoming neutrino results in a muon which points more 
closely to the initial neutrino direction at higher energies. The mean angular 
separation between an upward showering muon and its parent neutrino is 2.1° 
(assuming an atmospheric neutrino spectrum.) The corresponding numbers for 
upward through-going muons and upward stopping muons are approximately 
2.9° and 8.7° respectively. We also estimated the mean angular direction be- 
tween upward showering muons and their parent neutrinos for different parent 
neutrino energy spectra. For cosmic ray spectra and 1/E 2 spectra, the mean 
angular separation is 1.8° and 1.7° respectively. Therefore, we shall redo the 
searches for steady state point and diffuse astrophysical sources from Ref. [8] 
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with the showering muon dataset. There is no additional advantage in us- 
ing this dataset for transient astrophysical searches, since the timing window 
coincidence already provides a strong cut to reduce the background. 



8.1 WIMP searches 



We have performed searches for WIMP annihilations in the center of the Earth, 
Sun and Galactic Center using upward through-going muons [7]. Here we 
repeat the same search using upward showering muons. The cone size which 
contains most of the WIMP signal is inversely proportional to the WIMP 
mass. Since only high-mass WIMPs produce upward showering muons, we 
perform these searches in cones only up to 5°. Such a cone contains 90% of 
the signal for a WIMP of mass 1438 GeV from the Earth and 1000 GeV from 
the Sun and Galactic Center [7] . The observed data and expected background 
(evaluated in the same way as in Ref. [7]) are shown in Table 2. Since, there 
is no statistically significant excess in any of the search cones, we do not see 
any evidence for WIMP-induced upward showering muons in our dataset. 

Because of the reduced background from atmospheric neutrinos in this data 
sample, we expect to obtain better flux limits with this sample, compared 
to those in Ref. [7]. These flux limits calculated using only upward showering 
muons are plotted for the Earth, Sun and Galactic Center in Fig. 9, 10, and 11, 
respectively. The cutoff WIMP mass used for the calculation of WIMP flux lim- 
its using only the showering muon dataset is 1500 GeV for the Earth and the 
Galactic Center, and 2000 GeV for the Sun. This cutoff mass was calculated 
by determining the minimum mass for which the neutrino energy spectrum 
from WIMP annihilation is peaked at 1 TeV. This spectrum was calculated by 
assuming the same branching ratio in all the available annihilation channels 
for a given WIMP mass, and by using the analytic expressions for neutrino 
energy spectra calculated in Ref. [36] . The reason for the slightly higher cutoff 
for the Sun is because, unlike the Earth [37] and Galactic Center [38], ener- 
getic neutrinos lose energy or get absorbed by neutral and charged currents 
in the Sun [36,39,40,41,42]. To evaluate the WIMP flux limits with showering 
muons, we evaluated the showering efficiency at different WIMP masses. Note 
that we have also extended the flux limits with all upward muons down to 
10 GeV (from 18 GeV in Ref. [7]) by including the contributions of upward 
stopping muons. We have also shown for comparison the corresponding limits 
from AMANDA-II and Baikal for the Earth [43,44], AMANDA-II limits for 
the Sun [45], and MACRO limits for the Galactic Center [46]. The WIMP flux 
limits from all other detectors can be found in Ref. [7] . 
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Fig. 9. Super-K WIMP-induced flux limits of upward showering muons and all 
upward muons from the Earth as a function of mass. 
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Fig. 10. Super-K WIMP-induced flux limits of upward showering muons and all 
upward muons from the Sun as a function of mass. 

8.2 Search for a signal from potential point sources 



The equatorial coordinate distribution of upward showering muons is shown in 
Fig. 12. We look for signatures of neutrinos from 62 suspected point sources. 
Fifty-two of these sources were analyzed previously in Refs. [8,46] and many of 
these satisfy some of the features of the "beam dump model" [47,48] . Most of 
these sources are either supernova remnants, pulsars, magnetars or different 
types of active galactic nuclei. Some of these have been detected in TeV 7 
rays [49]. We have also considered some additional promising neutrino point 
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Fig. 11. Super-K WIMP-induced flux limits of upward showering muons and all 
upward muons from the Galactic Center as a function of mass. 



Table 2 

Observed and expected upward showering muons in cones with half-angles 3° and 
5° around the Earth, Sun and Galactic Center. 





Earth 


Sun 


Galactic Center 


Cone 


Data Background 


Data Background 


Data Background 


3° 


0.1 


0.2 


0.4 


5° 


0.5 


0.5 


0.8 



sources following the recent results from H.E.S.S 7-ray telescope [51,52]. 

In this search, we look for a statistically significant excess in a cone of angu- 
lar size 3° around these sources. Such a cone size would contain 90 % of the 
signal for a 1/E 2 neutrino spectra. The background was evaluated using at- 
mospheric neutrino Monte Carlo. Each Monte Carlo event was assigned a time 
sampled from the upward showering muon distribution, in order to match the 
livetime distribution of the observed events. The observed data and expected 
background are shown in Table 3. As we can see, there is no statistically signif- 
icant excess in any of these cones. For each source, we calculate the muon flux 
limit at 90% c.l. using the method described in Ref. [50]. To calculate the neu- 
trino flux limits, we assumed a E~ 2 neutrino energy spectrum. This spectrum 
arises in a number of astrophysical circumstances such as beam dumps and 
shock acceleration. The calculation of neutrino flux limits from the showering 
muon flux limits is as follows. 

The total flux of upward muons above an energy threshold E™ m for a 1/E 2 
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spectra is given by [47]: 

oo 

*v = A J dE u P,(E u ,E^)E~ 2 (9) 



/?min 



In Eqn. 9, A is a normalization constant for the differential neutrino flux, 
P M (E v , £™ in ) is the probability that a neutrino with energy E v produces 
a muon with energy greater than E^ m . The total flux of upward showering 



muons ($^ owerm 9) i s then given by 



= ^neW / dE™ ^~ mi ; ^ (iS™) (10) 



^ 



In Eqn. 10, e(P™ m ) is the probability for a muon to undergo radiative energy 
losses as a function of muon energy, k a i go is the efficiency of our algorithm to 
detect a true showering muon using the cuts in Sect. 5, and k cone is the fraction 
of the signal which falls within the 3° cone. For each source, we solve for A by 
substituting the obtained showering muon flux limits in the left-hand-side of 
Eqn. 10 and using $ M as evaluated from Eqn. 9. 



mm 
* ) 



To evaluate the above integrals, we calculated the values of P M [E v , E- 
using codes provided by M. Reno [29] (2005, private communication). More 
details on the assumptions used for this calculation can be found in Ref. [20]. 
The efficiency e(E™ m ) was calculated by first finding the fraction of muons 
from the mono-energetic Monte Carlo (used in Fig. 1) with dE/dX > 2.85 
MeV/cm at different muon energies, after which we applied a curve fit. The 
values of k a [ go and k cone were estimated for a 1/J? 2 spectra and each of them 
is equal to 0.9. Evaluating the integrals in Eqn. 10, the value of A for each 
source is given by A = 3.52 x 10 8 ^ owerin 9 GeV -1 . Once A is determined, 

oo 

the integrated neutrino flux limits were obtained by evaluating A J E~ 2 dE, 

where P™ m is the minimum neutrino energy which would make a showering 
muon and is chosen to be 10 GeV. This gives us 90% c.l. neutrino flux limits 
for a given source from the corresponding 90% c.l. showering muon flux limits. 
These muon and neutrino flux limits for all the sources are shown in Table 3. 
The corresponding neutrino flux limits for some of these sources from MACRO 
and AMANDA-II can be found in Refs [46] and Refs. [53,54]. A plot showing 
a comparison of the neutrino flux limits from these detectors for some sample 
sources as a function of declination can be found in Fig. 13. 
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Fig. 13. Neutrino flux limits from Super-K, AMANDA-II [54] (for E~ 2 neutrino 
spectra) and MACRO [46] (for E~ 2 neutrino spectra) as a function of declination 
for some selected sources. 

8.3 Search for unexpected point sources 



In order to search for a signal from an unexpected point source (which may not 
have any electromagnetic signature), we have looked for an excess of events 
within a cone of half- angle 4° around any upward showering muon event. The 
background is estimated using the atmospheric v Monte Carlo in the same 
way as in Sect. 6.1. The distribution of the number of observed events with 
the 4° cone is fit well by a Poisson distribution with a mean of 0.35 events. 
The comparison of data and Monte Carlo as well as results from the Poisson 
fit is shown in Fig. 14. Both methods show that the distribution of events is 
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Fig. 14. Number of events observed within a cone of half-angle 4° of any showering 
muon(dots). The solid line shows the expected atmospheric neutrino background. 
The dashed line shows a Poisson fit to the observed number of events with mean 
equal to 0.35 events 

consistent with the null hypothesis of having observed no signal from unknown 
point sources. 



8.4 Diffuse Searches from the Galactic Plane 

Although all the expected cosmic ray-induced neutrinos observed at Super-K 
arise from interactions in the Earth's atmosphere, some neutrinos could be 
produced from collisions of cosmic rays with hydrogen from the interstellar 
medium [55,56,57,58,59]. At low energies the expected flux from these neu- 
trinos is negligible compared to atmospheric neutrinos. According to some 
models [57], the minimum energy at which the atmospheric neutrino flux is 
equal to the flux of neutrinos produced from cosmic ray interactions in the 
ISM is about 10 TeV. Since most of the hydrogen is concentrated in the galac- 
tic plane, a signature of such events would be an excess from the galactic 
plane, which is characterized by galactic latitude of 0°. Searches for diffuse 
neutrinos from the galactic plane have been done with Super-K using all up- 
ward through-going muons [8] and with AMANDA-II [60]. We now repeat this 
search with the showering muon dataset. 

The galactic coordinate distribution of upward showering muons is shown in 
Fig. 15. As in Ref. [8], we looked for a statistically significant excess within 
±10° of the galactic plane. The observed signal events of 37 is consistent with 
the number of background events of 35, implying that there is no excess near 
galactic latitude of 0°. Thus we do not see any evidence that any of the upward 
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Fig. 15. Galactic latitude distribution of upward showering muons (dots) compared 
with Monte Carlo background from atmospheric neutrinos (solid line). 

showering muon events are coming from interaction products of cosmic rays 
with the interstellar medium. 
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Table 3. Table of observed number of upward showering muons and 
expected background from atmospheric neutrinos in a cone of half-angle 3° 
from selected point sources along with 90% c.l. muon and neutrino flux 
limits for a E~ 2 neutrino spectrum. 

Source Data Background // flux limits v flux limits 

(KT 15 cmr 2 s- 1 ) (llrtm-V 1 ) 
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Table 3 — Continued 



Source Data Background /i flux limits v flux limits 

(1Q- 15 cm- 2 s- 1 ) (lO^cm^s- 1 ) 
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Table 3 — Continued 



Source Data Background /j, flux limits v flux limits 

(lO-^cra- 2 ^- 1 ) (lO^cm-^- 1 ) 
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9 CONCLUSIONS 



From the sample of upward through-going muons in 1646 days of data we have 
isolated a sample of 318 showering muon events which lose energy through ra- 
diative processes. The mean parent neutrino energy of these events is ~ 1 TeV 
which make this subset the highest energy neutrinos seen in Super-K. The frac- 
tion of these events is approximately one fifth of upward through-going muon 
sample. At these high energies, the neutrino oscillation probability (for oscil- 
lation parameters obtained from Ref. [5]) is negligible. We have verified that 
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the zenith angle distribution of upward showering muons is consistent with 
null oscillations. This shows that the highest neutrino energy dataset (with 
mean parent neutrino energy of ~ 1 TeV, path-lengths of order 10000 km, 
and Am 2 ~ 0.0025 eV 2 ) is consistent with null oscillation. We also performed 
various searches for extra-terrestrial neutrino sources with this subset, such as 
from WIMP annihilations, suspected known and unknown point sources, as 
well as diffuse sources. We do not see any evidence for astrophysical neutrinos 
in the upward showering muon dataset. 
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